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Abstract 

We study the renormalization invariant trajectory of the <^> 4 -perturbation 
of the free field fixed point in the hierarchical approximation. We 
parametrized it by a running 4 -coupling g with linear step /3-function. 
We rigorously control the non-perturbative corrections to finite order 
approximants from double perturbation theory in g and g 2 ln(g). The 
construction uses a contraction mapping for the extended renormaliza- 
tion group composed of a hierarchical block spin transformation with 
a flow of g. 



1 Introduction 



The non-perturbative renormalization of c/> 4 -theory is a central problem of 
constructive quantum field theory p8|, [4C|, [7|]. The state of the art includes 



phase cell expansions [27, [L7], 35, |i[ |18| , renormalization group techniques 
g0|, |, |, ||, 0> H> fH IH> and random path representations §, |9|. How- 
ever, the non-perturbative renormalization of 4 -theory remains a mathe- 
matical enterprise of considerable difficulty. Recent work by Brydges, Di- 
mock, and Hurd § |, |l0 1 aims to simplify the older constructions and to 
cast renormalization theory into a more conceptual form. The present paper 
intends to make a modest contribution in the same direction. 

A key to the understanding of renormalization theory is Wilson's renor- 
malization group [45, 46 1 . The traditional starting point is a bare action 



So((f), 5o)Q- The goal is to compute a renormalized action as the limit n —> oo 
of S n ((j>,g n ) = R 7 l(So)(4>, go) of an iterated renormalization group transfor- 
mation Rl with scale L. The bare couplings go are tuned in this process so 
as to obtain a finite limit for g n . This process can be viewed as a trajectory 
in the dynamical system (on some space of actions) generated by Rl- If 
the action is a fixed point S*{(ft) of Rl then its renormalization becomes 
trivial: the bare action and the renormalized action become identical. This 
fixed point problem has a natural generalization. Consider a curve S((ft, g) 
of actions parametrized by a (running) coupling g such that 

1. S(<p, 0) is a fixed point of Rl, 

2. d g S(4>, g)\ g =o is an eigenvector O(cft) of the linearization of Rl at S+((ft), 
and 

3. Rl(S)(((), g) = S((f>,5^ 1 (g)), where 6l is a step /^-function. 

In this case, the bare action and the renormalized action have the same 
functional dependence of (ft but correspond to different values of the running 
coupling g. Renormalization amounts to control the flow generated by the 
step /3-function, a comparatively easy task. 

We will study a variant of this renormalization problem, where Rl is a 
block spin transformation for a three dimensional scalar lattice field theory 
with hierarchical covariance [32]. It is designed such that the interaction 
remains local under the renormalization group evolution. The lattice inter- 
action Boltzmann factor factorizes into product of local Boltzmann factors 
Z(<ft) = exp(— V{4>j) (one for each lattice site), which are functions of a real 
variable (ft. In three dimensions, the hierarchical renormalization group then 
reduces to the non-linear integral transformation 

R L (z)W = jy dMO z (V^ + c) } , (i) 



lr The quantum field is understood to be rescaled to a unit ultraviolet cutoff. 



1 



where d//(£) is the Gaussian measure on R with mean zero and unit covari- 
ance. This transformation and variants of it have been studied by many 
authors both as a model of constructive renormalization and also because of 
its properties as a non-linear theory. Rigorous work on hierarchical models 
(more generally ultralocal renormalization groups) includes: 

1. the e-expansion [^] 



2. the 03 infrared fixed point |[l|, 33, 34 

3. the massive perturbation of the 0|. fixed point [32] 

4. the renormalization group differential equation 

5. the 04 infrared problem |37], [| 

6. the 04 ultraviolet problem at negative coupling [24] 

7. the renormalized (/^-trajectory p3] 



8. the S£/(2)-lattice gauge theory |3| 

9. the non- linear a- model |26|, [38| 

10. the sine-Gordon model []l4|, 3C] 

11. multigrid expansions |37], |3^[ 

12. and random surfaces ll 



Beyond the hierarchical approximation, one has to deal with non-local in- 
teractions generated by the renormalization group. Although non-local cor- 
rections are rather small in all models brought under control so far, the 
mathematical apparatus needed to control them is a lot more sophisticated, 
the main tool being polymer expansions. The virtue of hierarchical models 
is that they allow to study renormalization effects without this additional 
burden (or perhaps joy). 

In this paper, we continue the work started in |44]] . We look for a curve 
of renormalized interaction Boltzmann factors Z((f),g) with the following 
properties: 

1. Z(4>,g) = exp(— g : 4 :) ^1 + 0(g 2 ^j , i.e., Z((f>,g) emerges from 

the trivial fixed point Z*(0) = 1 (the free massless hierarchical field) 
tangent to a (normal ordered) 4 -vertex, and 

2. Rl(Z)(i/j, 5l{9)) = Z(ip,g), i.e., Z((f>,g) is a fixed point of the ex- 
tended renormalization group Rl X 5^ with a linear step /^-function 

S(g)=L- 1 g, 
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The problem is thus to construct a non-trivial fixed point of the extended 
renormalization group Sl = Rl x $l- We do this by means of a contraction 
mapping. For this purpose, we split Z = Zi+Z 2 , where Z\ is an approximate 
fixed point, and where Z 2 is a correction. We iterate the transformation of Z 2 
with Z\ kept fixed. This transformation is shown to contract, provided that 
Z\ is in a certain sense a sufficiently good approximation. We compute Z\ as 
a polynomial approximant of finite order in a (formal) double perturbation 
expansion in g and g 2 ln(<?). We then prove that this approximant is indeed 
sufficiently accurate provided that the order of perturbation theory is at 
least seven. The result of this construction is the following Theorem.^ 



Theorem 1.1 Let F 2 {g) be a continuous positive function of the form 



F 2 (g) 



9 exp(c*0' 



+ cg) 



(2) 



with positive constants a, c+, and c. Let ZQjj((p,g) be the quadratic fixed 
point 



Zqu(4>,9) = exp a QU (g) 



bQu(g) 



of Sl. Let B 2 be the Banach space of functions Z 2 ■ 
respect to the norm 



x [0, g n 



\Z 2 \\f 2 



sup 

(<t>,g)€Rx[0,g n 



ZQu{(t>,g)F 2 {g) 



(3) 



with 



(4) 



LetZi((j),g) = exp (— W rmax )(c/>, g)) , where V^ r ' m ^ '(</> , g) is the polynomial 
approximant of order r max (in g) of the perturbative solution to the fixed 
point equation as a double expansion in g and g 2 ln(g) ■ 



For r„ 



7, there exist positive constants g max , a, 



and Cr> such 



that the transformation Sl {Z\ 
ball 



Z 2 ) — Z\ is a contraction mapping on the 



{Z 2 €B 2 \ \\Z 2 \\ F2 <C 2 }. 



(5) 



It follows that there exists a unique fixed point in this ball. Furthermore, 
the iteration of the contraction mapping gives a convergent representation 
for this fixed point. 

The renormalized ^-trajectory was constructed in for all dimen- 
sions 2 < D < 4 with the exception of a discrete set of special dimensions, 
where resonances of power counting factors occur ]4"H| . Unfortunately, the 
case D = 3 is such a resonant case, and was therefore excluded in [441. The 



problem is that our renormalization problem does not have a formal power 



2 Zqu and V"' rmax ' are explained in the bulk of this paper. 
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series solution in g in three dimensions. However, it does have a solution as 
a formal double perturbation expansion in g and g 2 ln(g). The main con- 
tent of this paper is to deal with these logarithmic corrections. Polynomial 
approximants from this double perturbation theory turn out to suffice for 
the contraction mapping. The backbone of our approach is the contraction 
mapping. This part is identical in resonant and non-resonant dimensions. 
To keep this paper selfcontained we have included a section on the contrac- 
tion mapping with fresh proofs and improved bounds as compared to 44]. 
In particular, we present 

1. a better scheme independent proof of the contraction property, 

2. an example of bounds, which are true for all couplings and not only 
small couplings, 

3. a better and more explicit treatment of the tree approximation, 

4. and last not least a full stability analysis of the g-g 2 ln(<?)-approximants. 



Unlike [44], analyticity in cj) is not used here. This paper is organized as 
follows. Section two contains a brief review of the hierarchical renormaliza- 
tion group. Section three is devoted to the contraction mapping method. In 
Section four, we prove a stability bound and an error bound for the first or- 
der approximant. It serves as a template for the higher order approximants, 
which are analyzed in Section five. We conclude with a few remarks and 
outlooks. 



2 Hierarchical renormalization group 

The hierarchical renormalization group is the theory of a non-linear integral 
transformation R acting on a certain space of functions Z : M — > R. We 
speak of Z as a hierarchical model in Statistical Physics if Z is positive, 
continuous, and in a certain sense measurable. 



2.1 Renormalization group transformation 

Hierarchical renormalization groups come in different versions. We mention 
the transformations of Dyson [15], Wilson [45], Baker |2|,and Gallavotti 
p0| . The latter version has been investigated as a model for the constructive 
renormalization of massless scalar field theories by Gawedzki and Kupiainen 

] . Following 



|23], [2J, by Pordt |37|, gg], and by Koch and Wittwer |31 
their path, we define R as follows. 

Definition 2.1 Let L denote a positive integer, L £ {2,3,4,...}, and let 
D = 3. Let a = L D , (3 = L 1_ ~, and 7 = 1. Let R be the non-linear 
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integral transformation given by 



R{Z){^) = I J d^cozc^ + oj , (6) 

where d// 7 (£) denotes the Gaussian measure^ on R urai/i mean zero and co- 
variance 7. 

L is the block scale. Eventually, we will choose L to be large. D is the 
dimension of the (Euclidean) space-time. We will restrict our attention to 
the three dimensional transformation. In the sequel, D and 7 are constant, 
D = 3 and 7 = 1. The transformation (g) calls to be supplemented by a 
domain. The following statement defines our model space. 

Definition 2.2 Let M be the space of continuous even functions Z : R — ► R 
with finite norm ||Z[|oo = su P</>ell |Z($)| • 

A part of our analysis will be to restrict -R to suitable invariant subspaces 
of M, for instance subspaces of functions with a rapid decrease at infinity. 
Positive functions form a subset A4 + of Ai. But 7W + is not a linear subspace 
of M. Because we intend to use Banach space theory, we take the latter as 
our starting point. 

Lemma 2.1 The transformation R acts on the space M. It satisfies the 
bound H-RCZJHoc < \\Z\\% . 

2.2 Trivial fixed point and tangent map 

R has two trivial non-zero fixed points in A4 + , the ultraviolet fixed point 
Zuv((p) = 1 j an d the (Gaussian) high temperature fixed point 

Z HT (cl>) = AHre-^* 2 , (7) 

where 



afi 2 



A H t = (a/3^) 2 ^- 1 ) = Li 7 ^ 11 , b HT = — = L z - 1. (8) 

7 

In this paper, we will study the unstable manifold of the ultraviolet fixed 
point. Basic knowledge about it is gained from the linearized renormaliza- 
tion group. Recall the following well known facts. 

Lemma 2.2 Let O G R[<j) 2 ] be an even polynomial in cj) with coefficients in 
R of the form 0(<p) = 4> 2n + lower powers of<p 2 . Let Z : R x R+ -> R be 
defined by 

Z(4>,g) = e-»°W. (9) 

3 The Fourier transform of the Gaussian measure is f dp 7 (C) e 1 ^ = e -7 * 1/2 . 
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(I) For all g G R + , Z(-,g) £ A4+ ; and Z(</),0) = Z uv (4>). (II) For all g £ 
M + , R(Z)(-,g) is continuously differentiable in g, and lim g ^Q+ d g R(Z)(-,g) 
defines a linear operator 

D Zuv (R)(0)(4>) = a [ d^ 7 (C) <W + 0- (10) 



on 



The linear operator Dz uv (R) is the tangent map of R at Zjjy This tangent 
map is diagonalizable. 



Lemma 2.3 Let H2n be the Hermite polynomial of order 2n j^9\ 8.95]. Let 
P2n £ M[</> 2 ] be given by 

«*<*.») = (f)"H 2 „(^=), " = ^- (ID 

M^e /iai>e £/ia£ Dz uv (R)(P2n) = ^2n.P2n, where \2 n = cef3 2n = L U2n with 
scaling dimension a<in = D + n{2 — D) = 3 — n. 

In three dimensions, Pq, P2, and -P4 are relevant, Pg is marginal, and all 
others are irrelevant. From the set of eigenfunctions, we select the relevant 
non-quadratic member P4. Its explicit form is 

P A (<f>,v) = cf) A -6v(j) 2 + 3v 2 , cj 4 = A-D = 1. (12) 

Another way of saying that P4 is an eigenfunction of the tangent map with 
eigenvalue A4 is the following. Let V(4>, g) = g P^((p,v). Then we have that 

D Zuv (R)(V)^,Sg) = Vty,g), 5' 1 = af3 4 = L^. (13) 

In other words, V((p,g) is a fixed point of Dz uv (R) X 5*, the tangent map 
extended by a flow of g. The linear function g — > 5g will serve as our step 
/3-function. 

As in [44], the goal of this paper is to construct a non- linear analogue of 



(13), which is a restricted homeomorphism from the tangent space at Z\jy 
to an invariant cone in A4, originating at Zjjv- 



2.3 Extended renormalization group 

From now on, we turn our attention from points in A4 to parametrized 
curves in M., originating at Z\jy. 

Definition 2.3 Let <? ma x £ ^ + and G = [0,g max ]- Let N be the space of 
continuous functions Z:txG->K such that Z{(p,g) = Z(—(j),g) for all 
{4>,g) ElxG, Z(4>,0) = Z uv (4>), and \\Z\loa = sup {MmxG \Z(4>, g)\ is 
finite. 
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For technical reasons, we introduced a maximal coupling g max - By default, 
this maximal coupling is an arbitrary large number in the following. Re- 
strictions on g max will be explicitely stated. 

Definition 2.4 Let 04 = 4 — D and 5 = L~°~ 4 . Let S be the non-linear 
transformation given by 

S(Z)ty,g) = (Rx5*)(Z)(iP,g) = R(Z)ty,6g). (14) 



Lemma 2.4 The transformation S acts on the space M. It satisfies the 
bound IIS^Hoo < \\Z\\^. 

This bound is an immediate consequence of SG C G (since Sg < g) and 



sup \S(Z)(<p,g)\ < { sup \Z(<f>,g)\\ (15) 
(<t>,g)eRxG [(<t>,g)mx8& J 

Functions in M are bounded but not necessarily decaying at large fields. To 
bound non-perturbative contributions, we will need an exponential decay 
at large fields. The following observation suggests how to select a subspace 
Mqu C J\f, which suits this purpose. 

Lemma 2.5 Let Zqjj : R x G — > M. + be the positive continuous function 
given by 

Z QU (cP,g) = exp \a QU {g) - ^|^0 2 | , (16) 

where 

a Q u(9) = ^ £« W1 °g{ i + g / \ ^ 

n=l 

and 

a R 2 - 1 q p 

bQU (g) = 2E. (18) 

where p = = 2. The function Zqu is an element of M. It is a fixed 

point of S. 

The assignment g 1— ► ZQu(-,g) is a continuous parametrized curve in 
There is no reason, not to set G = R + at this point (or g max = 00). Then 
Zqu becomes a curve of Gauss-functions, which connects the two trivial 
fixed points of R, the ultraviolet fixed point Zqu((/>,0) = Z\jy{§) to the 
high temperature fixed point Zqu(4>,oo) = Zht{4>)- 
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Definition 2.5 Let Nqu be the subspace of N , consisting of functions with 
finite norm 



\Z\ 



QU 



sup 

(0,p)eKxG 



z Qu(4>,g) 



(19) 



completed to a Banach space. 



Functions in Nqu are in particular continuous in g. Their most useful 
property is the inbuilt bound 



\Z( 



< \\z\ 



QU ZQU 



(20) 



At fixed g, it compares the decay of Z at large <f> with the decay of the fixed 
point Zqjj. (^) serves as the basic large field bound in this paper. 

Lemma 2.6 The non-linear transformation S acts on Nqu- (It leaves in- 
variant Nqu C N.) Let Ze Nqu- Then \\S{Z)\\ QU < {{Z^. 

We can now state our goal more precisely as to find a non-trivial fixed point 
of S in Nqu (besides the ^-independent trivial fixed points of R, and besides 
Zqjj). To this aim, we need to introduce yet another subspace oINqu- 

Definition 2.6 Let g max < oo. Let F : G — > R + be a continuous positive 
function with finite norm ||-F||oo = su P 9 eG l-^XsOl • Let Nf be the subspace 
of Nqu consisting of functions with finite norm 



\Z\ 



sup 

0,g)eRxG 



Z(<P,g) 



Z QU (ct>,g)F(g) 



(21) 



(|2lf) refines (|20|). Its motive is to gain control over the ^-dependence of 
ZJQ It will be crucial to carefully choose the function F. To give a flavour 
of, what kind of function F should be, notice the following fact. Z G Nf 
implies that 



and thus 



\Z(<l>,g)\ < \\Z\\ F Z QU (<t>,g) F(g). 



\S(Z)(<j>,g)\ < \\Z\\ F Z Q u(<j ) ,g)F(5gy 



(22) 



(23) 



If F is such that the condition F{5 g) a < F(g) holds for all j£G, then Nf 
is invariant under S. Furthermore, the ball H^Hf < 1 is then mapped to 
itself under S. Perturbation theory suggests functions of the form F(g) = 
C g a exp(cg T ). Depending on the actual values of C, a, and r, the above 
condition may pose a restriction on <7 ma x- 

4 The information that Z is an element of Nqu implies no more information about its 
^-dependence than that, for all <j> € R, Z(cj>, g) is a bounded function of g. 
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3 Contraction mapping 



In this section, we present our method to construct a non-trivial fixed point 



of S. As in 1 44], we split Z = Z\ + Z 2 into two parts Z\ and Z 2 , where Z\ 
is an approximate fixed point, and where Z 2 is an error term. Z\ will be 
kept fixed. We intend to estimate the non-linear mapping of the Z 2 under 
certain assumptions on Z\. 

Lemma 3.1 Let Z\ and Z 2 be two elements of Nqij. Then 

S(Z 1 + Z 2 )-Z 1 = T 1 (Z 1 )+T 2 (Z 1 ,Z 2 ), (24) 



where 



and 



Ti(Zi) = S(Z 1 )-Z 1 (25) 



T 2 (Z 1 ,Z 2 ) = S(Z 1 + Z 2 ) - 5(Zi). (26) 

Since S acts on Mqu, the mappings T\ and T 2 are both well defined. The 
norm of Ti(Z\) measures the quality of the approximate fixed point Z\. 
Equipped with a bound on Ti(Zi), the issue is to find bounds on Z 2 which 
iterate under (U). 

3.1 Invariant ball B 

The dependence of Z\ and Z 2 on g will be controlled with two different 
functions F\ and F 2 to be specified below. 

Lemma 3.2 Let Fi : G — > M + , i = 1,2, be two continuous positive func- 
tions. Abbreviate Mi = A/"^ and \\Zi\\ = \\Zi\\p.. Assume two functions Z^ 
with Zi &Ni. Thus for all ((f), g) £ R x G, 

1^(0,5)1 < \\Z i \\Z QU {<l>,g)F i {g). (27) 

Then T 2 {Z\, Z 2 ) is an element of Mqu- For all (4>,g) eRxG we have that 

\T 2 (Z x ,Z 2 \4>,g)\ < 

\(\\Z 1 \\F 1 {6g) + \\Z2\\F 2 (Sg)) - ( ||Zi|| F^Sg) ) ) Z QU (<f>,g). (28) 
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Proof T 2 (Zx, Z 2 ) has the integral representation 



T 2 (Z 1 ,Z 2 )(^,g) 



J dt-S(Z 1 +tZ 2 )^,g) 
J^dt^lJ d^ 7 (C) (Z 1 + tZ 2 )(P^ + C,6g)\ 
X dta { I d H {Q (Zx + tZ^i^ + CSg)^ 



a-l 



x J d^(0 W + C,«ff) (29) 

Using the bound (p7|), it follows that 
\T 2 {Z x ,Z 2 )(^,g)\ 

< J dt a (\\Z x \\F 1 {8g)+t \\Z 2 \\ F 2 (6g)J 

x \\Z 2 \\ F 2 (5g) U d^ 7 (C) Z QU {M + t,5g) 

= J Q ^^{\\Zi\\F x [8g)+t\\Z 2 \\F % {Sg)\ Z QU {^,g) 

= y\\Z 1 \\F 1 (5g) + \\Z 2 \\F 2 (5g)^ - F 1 (6gfj | 

x Z QU (,/>,g). □ (30) 



Differing from |23[ |, no split into small and large fields is required here. 
Small and large fields are taken care of simultaneously by the ^-dependence 
of Zqxj. For the right hand side of (p8|), we notice the elementary estimate 



a-l 



a [|Zi|| FxiSg) + \\Z 2 \\ F 2 (5g) \\Z 2 \\ F 2 (5g). (31) 



Lemma 3.3 Let Z\ be as in Lemma 3.L Let Z\ be such that in addition 



||Ti(Zi)|| = ||Ti(Zi)||i? 2 is finite. Then we have that for all ((f), g) ElxG, 

1^(^X0,3)1 < \\T 1 (Z 1 )\\Z QU {4>,g)F 2 {g). (32) 

Let n be an integer, n £ {1, 2, 3, . . . }. Assume that F\ and F 2 conspire such 
that 

a(\\Z 1 \\F 1 {8g) + {n+l)\\T l {Z l )\\F 2 (5g)\ F 2 (5g) < F 2 (g). 
V / n+l 

(33) 
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Let B be the ball in M 2 given by 

\\Z 2 \\ < (n + 1) \\mZt) 



(34) 



Then it follows that B is invariant under the transformation (24)- For all 
Z 2 £M and for all (</>, g) G R x G, 



\S(Z 1 +Z 2 )(^g)-Z 1 (<j ) ,g)\ < (n+1) {{T^W Z QU (<f>,g) F 2 (<f>,g). 



(35) 



Proof From (||), (HJ), (||), and (||), it follows that 

\\S(Z 1 + Z 2 )-Z 1 \\ < ||Ti(Z 1 )|| + ||T 2 (Zi,Z 2 )|| 
< [|Ti(Z0||+ " 



n+1 



Zo 



< \\Ti(Z 1 )\\ + 



n 



n+1 

(n + l)||Ti(Zi)||. □ 



(n + 1) wmzt 



(36) 



We remark that the additional assumptions made in Lemma |3.3| do not clash 



with those made in Lemma 3.2 



Lemma 3.3 suggests the following strategy. We look for functions Z\, F\, 
and F 2 such that H^iH^ and ||Ti(Zi)||^ 2 are both finite. Furthermore, F\ 
and F 2 have to satisfy (^) . Then we have an invariant ball of error terms 
Z 2 in the norm associated with F 2 . 

3.2 Contraction mapping 

Our next task is to establish the contraction property for the transformation 
(H) on the ball B. 



Lemma 3.4 Let Z\, Z 2 , and Z' 2 be as in Lemma 3.L Then we have for all 
(</>, g) G R X Q the bound 



\S{Z l + Z' 2 ){cj>,g)-S{Z 1 + Z 2 ){<j>,g)\ < 

a(\\Z 1 \\F 1 {5g)+ sup \\Z 2 + t (Z' 2 - Z 2 )\\ F 2 (5g)) 
V te[o,i] / 

x Z QU (ip,g) \\Z' 2 - Z 2 \\ 



a-l 



F 2 (5g) 



(37) 
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Proof From the integral representation 

S(Z l + Z , 2 )(1>,g)-S(Z 1 + Z 2 )y>,g) 
r 1 £i 

= J dt-S(Z 1 + Z 2 + t (Z' 2 - Z 2 )) (V, g) 

= J dt §- t {J d ^(0 ( z ± + Z 2 + t{Z' 2 - Z 2 )) + C, Sg) 

= JdtaU d^ 7 (C) [Z x + Z 2 + t{Z' 2 - Z 2 )) (04, + C, Sg) 

x j d^ 7 (C) (Z' 2 -Z 2 )(J3if, + C,6g) 

it follows that 

\S(Z 1 + Z' 2 )(tp,g) - S(Z 1 + Z 2 ){il>,g)\ 



a-l 



(38) 



< J dta \\\Z 1 \\F 1 (Sg) + \\Z 2 + t{Z' 2 - Z 2 )\\ F 2 (5g) 
x F 2 {5g) { [ o> 7 (C) Z QU (W + (,5g)\ \\Z' 2 - Z 2 \ 



a-l 



<a(\\Z 1 \\F 1 (5g) + sup \\Z 2 + t{Z' 2 -Z 2 )\\ F 2 (5g)) F 2 (Sg) 

V t€[0,l] J 

x Z QU (iP,g) \\Z' 2 -Z 2 \\. □ (39) 
The contraction property follows from this estimate in cooperation with 



Lemma 3.3 



Lemma 3.5 Let Z\, Z 2 , and Z' 2 be as in Lemma \3. % For all Z 2 6 B ana 

Z 2 £ B, we have that 



\\S{Z l +Z' 2 )-S{Z 1 + Z 2 )\\ < 



n 



n + l 



\z 2 ~ ZvW- 



(40) 



Proof Since B is convex, it follows that for all t £ [0, 1], 

\\Z 2 + t{Z' 2 - Z 2 )\\ F 2 (5g) < (n + l) ||Ti(Zi)||. (41) 
We can therefore use ( |3"4| ) to conclude the asserted bound from (|37|) , namely 
\S(Z 1 + Z' 2 )(4>,g) - S(Z 1 + Z 2 )(<P,g)\ 

71 

< —— Z QU {4>,g) F 2 (g) \\Z> 2 - Z 2 \\. □ (42) 
n+l 



Corollary 3.1 Let B be as in Lemma 3.S. B is invariant under the trans 



formation (Ww- The transformation n2A) is a contraction mapping. The 



transformation (2£) has a unique fixed point in B. 
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The iteration of the contraction mapping is a convergent scheme to compute 
the fixed point. Let Z2.0 = and define a sequence of (non-perturbative) 
approximants by 

Z 2 , n+1 = S(Zi + Z 2>n ) - Z x . (43) 

As n — ► 00, this sequence converges to the desired fixed point. We cannot 
say much about this limit at this level of generality. However, if Z\ is strictly 
positive than it follows immediately that the fixed point is non-negative. 



3.3 Estimates for all couplings 

Our next subject is to construct examples of the two functions F\ and F 2 
with the property (33). These examples are tailor made for perturbation 
theory. 

Suppose that we are told an approximate fixed point Z\ from an inde- 
pendent calculation, whose large field decay is (at least) Gaussian. Then we 
proceed as follows. 



1. We determine F\ such that 



sup 



Zi((f),g) 



ZQu(<t>,g) 



< F 1 (g). 



(44) 



Then ||^x||fi < 1- Since it contains information about the large field 
behavior of Z\, we call this bound the stability bound on Z\. 



2. We determine F2 from an estimate on Ti(Zi) such that 

Ti{Z x ){4>,g) 



sup 



ZQu(4>,g) 



< Ci F 2 (g) 



(45) 



for some finite constant C\. Then ||Ti(Zi)||i? 2 < C\. We call this 
second bound the error bound on Z\. 

3. We check that F\ and F 2 satisfy (]3 

If F\ and F 2 are functions of the type of this section, then the validity of 
( |33D is guaranteed by the below estimates. 

Remember that Z\ plays the role of an approximate fixed point. This 
means that T\{Z\) should be small compared with Z\. A construction for 
all couplings should in particular cover the case of small couplings. Inspired 
by perturbation theory, we choose 



F 2 {g) = g a F x (g). 



(46) 



It says that the error term goes to zero as the coupling goes to zero. The 
speed of this process is given by the exponent a. 
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Lemma 3.6 Let F\ : G — > M + be a positive continuous function. Let 
II^iIIfi < 1- Let a be a positve real number, with a > a*, wher^ 

°* = 4^ = 3- (47) 

Let F 2 (g) = g a F\{g)- Let \\Ti(Z{)\\f 2 < C\ for some positive constant C\. 
Let n £ N. Put C 2 = (n+ 1)C±. Assume that F\ has the following property. 
Let there exists a positive constant Cp such that for all g £ G, 

e ( a -l)C 2 (Sgr Ft(Sg) a < C F F 1 {g). (48) 

Then we have the following estimate. There exists a positive real number 
L m ; n such that for L larger than this number, we have the bound 

a i\\Z 1 \\F 1 (8 g) + (n+ I) \\T 1 {Z 1 )\\ F 2 [5 g) J F 2 (S g) < — [ F 2 (g). 

(49) 

Proof From the assumptions (|45|), fl46|), and (fi"8|), it follows that 
a (\\Z 1 \\F 1 (5g) + (n + l)\\T 1 (Z 1 )\\F 2 (5g)] F 2 {5g) 



< ol (l + C 2 (Sg)A {Sgf F^Sgf 

< a^ a - 1)C2 ^ r (Sgf F^Sgf 

< a (5gf C F F 1 (g) 

= ad* C F F 2 (g). (50) 

For a > a*, the L-dependent factor aS a = L D ^ D ~ can be made 
arbitrary small by taking L to be large. □ 

We learn that ( p3| ) (identical with fl49|)) holds provided that F\ satisfies (48). 
An example of a function F\ which satisfies (E3) without a restriction on 
<?max is the following. 

Lemma 3.7 Let c* and c be positive constants. Let a be a positive constant 
such that a > a*. Let F\ be the function 

Ft(g) = exp(c*f* + cg a ). (51) 
For all positive constants C 2 such that 



5 The exponent cr* is such that aS a * = 1. For a > a+ the block volume a is beaten by 
the power of the step /3-function. 
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this function Fi satisfies the bound 

e ( a -l)C 2 (5g)° F ^ S g) a < Fx(g). (53) 

Proof For this particular function F\, we have that 

exp{(a-l)C 2 (S g y} exp{c* 5 CT * + ac(<5< ? )< 7 } 

< enp{c*g a * + a5 a (C 2 + c)g (T } 

< exp{c*g a * + cg a }. □ (54) 



Lemma 3.7 contains a restriction on C%. To avoid a conflict with fl33|), the 
constant c has to be chosen to be sufficiently large. 

The pair Fi given by ( |5l| ) and (^) (with a > <j*) satisfies all properties 
needed for the contraction mapping. Remarkably, it involves no restriction 
on the size of <7 max - With this pair, we could construct the </> 4 -theory at 
arbitrary large couplings. The problem is to compute an approximate fixed 
point Z\, which satisfies the stability bound (f44|) with this function F\, and 
which satisfies the error bound with a sufficiently large exponent a. Un- 
fortunately, we have not succeeded to prove these bounds for approximants 
from perturbation theory. 

3.4 Estimates for small couplings 

Therefore, we supplement the bounds for all couplings by suitable bounds 
for small couplings. There are two sources of constraints on the value of 
Smax- One source is that we may not be able to prove the stability bound on 
Z\ for arbitrary large couplings. Also we may not be able to prove the error 
bound on T\(Zi) for arbitrary large couplings. In this section, we will not 
speak about this source of problems. These constraints cannot be addressed 
before we actually compute Z\. The second source is that the function Fi, 
which we find from ([5]) (say by defining F\ by equality) , might not meet the 
requirement ( |48| ) for arbitrary large couplings. If both effects come together, 
we have to put <7 max equal to the minimum from these constraints. 

Lemma 3.8 Let c* and c be a positive constants. Let t be a positive con- 
stant such that t < a*. Let F\ be the function 

F 1 (g) = exp( C * 5 CT * +cg T ). (55) 

Let C2 andCF be positive constants, withCp > 1. For all values of L , there 
exists a maximal coupling g(L) such that, for all g less than this maximal 
coupling, we have that 

e ( a -l)C 2 (5 9 r F x (8g) a < C F Fx(g). (56) 
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Proof The estimate is very similar to the one for all couplings. The identity 



e (a-l)C 2 (6 9 y Fl ^g) a = e <=(a* T -i)<r F^g) (57) 
yields the condition 

e (a-l)C2(5g) a +c(a5 r -l)g- r < Q p 

on g. Given exponents a and r, together with constants C2 and c, and a 
scale L, (]58|) determines g(L). □ 

To get an idea how g(L) behaves as a function of L, we neglect the the first 
term in (|58|). Then 

9(£) < ( J^U } (59) 



c (a5 T - 1) 

shows that g(L) shrinks as a certain power of L. A typical value of r in 
our perturbation theory is one. Since <r* = 3 in three dimensions, we are 
in the small coupling case. The good news is that the maximal coupling is 
not necessarily ridiculously small. 

It is likely that there exist other functions Fi which meet the require- 
ments of the above contraction mapping. The particular ones considered 
here are made for approximants Z\ from perturbation theory. The inter- 
esting question remains whether this contraction mapping works with other 
approximants, for instance approximants from numerical work. This would 
be more in the spirit of the fixed point construction of [31]. 



4 Linear approximation 

The simplest approximation is a pure c/> 4 -vertex. In this section, we will 
prove a stability bound and an error bound for this linear approximation 
defined by 

Z(c/>,g) = e- v ^\ V&g) = gP^,v). (60) 

(To simplify the notation, we write Z instead of Z\.) The linear approxi- 
mation will not suffice to obtain a contraction mapping in three dimensions. 
But it is a part of, and also an instructive example for, the finer bounds to 
be presented below. 



4.1 Stability bound 

The classical stability bound for (|6^) relies on analyticity in the field variable 
a strip around the real axis [^3| . We proceed differently therefrom |37], . 
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Lemma 4.1 The polynomial P^<f),v) is bounded from below by 

P^,v) > ^ 4 -15zA 



(61) 



Proof 



b 4 - 6vcj) 2 + 3v 2 
^ 2 -^) 2 -(^) 2 + (l- e 2 )0 4 + 3. 2 



> (1 - e 2 ) <f + 3 ( 1 
For e 2 = i, the assertion follows. □ 



(62) 



Lemma 4.2 For all (4>,g) GlxG, we have that ft6(\ ) is bounded from above 
by 



Z(<P,g) < Z QU (0,g)e a ^ 



(63) 



with 



ai{g) < 15v 2 g + lg 2p ~ 1 -a QU (g). 



Proof For any c, we have the elementary bound 



(64) 



> 4 = ^"f 



)2 + c0 2_^ > Cf ,2 

4 



Put c = s^ -1 to obt 



am 



(65) 



g^ > 9 P A2 9 2p ~ x 



2 T ~ 2 T 8 
The values of L, D, and 7 are such that ( |i~8l) is bounded from above by 



(66) 



&Qt/(s) = rn.. , - n < — 9 P < 9 P - □ 



L D 7 1 + 5^ 



7 



(67) 



The bound (p3) suggests a function Fi of the form (51). 
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Lemma 4.3 Let Z be given by ([figy. Let F\ : G — > M + be the function 

= exp 15 v 2 g + ?-- 



(68) 



Then Z is bounded in the norm (21). We have that \\Z\\p 1 < 1- 
Proof For all g > 0, aQjj(g) > 0, wherefore 



l-Fl 



sup 

(0,g)eKxG 



ZQu((f>,g)Fi(g) 



< sup e ~ a< 3 u(9) < 1. □ (69) 



In three dimensions, the value of p is two and that of <r* is three. (Recall 
(18) and (^).) Coincidentally, 2p — 1 = er*. But the second term in the 
exponent in ( |68| ) is only linear in g. Although the stability bound ( |68| ) holds 
for any value of g, it restricts our contraction mapping to small couplings. 



4.2 Error bound 

Equipped with this stability bound, one is led to estimate T\(Z) in the norm 



given by (46), where a is a suitable exponent. As we will see, in fact as we 
know from 44], such a bound indeed holds. But the exponent a is only 
one half and therefore smaller than a*. The linear approximation therefore 
suffices only for a construction in low dimensions. Let us nevertheless see 
how the exponent one half comes about. For this purpose, we consider the 
following interpolation formula. 

Definition 4.1 Let X : R x G x [0, 1] -> R be defined by 

Xy,,g,t) = lj d Mt7 (C) exp(- J d/i (1 _ i)7 (£) H^ + C + j 

(70) 

In the limit of a vanishing covariance, the Gaussian measure d// 7 (£) be- 
comes a Dirac measure d£<5(£). Therefore, (|69|) interpolates between the 
exponentiated linerarized renormalization group transformation ([lC|) 



X(tf,g,0) = exp(^-« j (^(OW + f.iff) 



and the full renormalization group transformation 

X&,g,l) = { [ dny(C) exp(-V(W + (,6g) 



(71) 



(72) 



both transformations being extended by the flow of g. The usefulness of this 
interpolation relies on the following property. 
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Lemma 4.4 For t £ (0, 1), X is continuously differentiable in t. We have 
that 





—X(^,g,t) = 
a I J d^ 7 (C) exp (- J dM(i-t) 7 (0 *W + C + e, <5 5 ) 
X |y d/x t7 (C) exp ^- y d/i (1 _ t) 7 (f ) + C + C «5) 

x 2 (I; / d ^d-*) 7(0 *W + c + e, <55) 



a-1 



(73) 



Proof 



^ y d Mi7 (o expf- y dM(i-t )7 (o n^+c+o) 

"7 <9 2 9 



d/At 7 (C) 



2 3C 2 + 0* 



7 d 2 9 
2 5C 2 + St. 



(74) 



and 



7 3 2 9 

2 sc 2 + at 



d/x (1 _t )7 (e)F(^ + C + C) = 0. 



□ 



(75) 



Notice that (73) is of second order in V due to the cancellation (fT5|) . Notice 



furthermore that fl73| ) is non- negative. The integral of (73) yields a represen- 
tation for T\(Z), which can be used to derive an upper bound of the desired 
form. 



Lemma 4.5 Let Z be given by \6(\). Let X be given by (uQ). Then we have 
that 



/*1 ^ 

T x (Z){^g) = X(i>,g,l)-X(i>,g,0) = J dt — Xii.'.g. I). (76) 
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Proof Because V is an eigenvector of the extended linearized renormaliza- 
tion group (|T^), 



a J W^ + f.is) = V(il>,g) (77) 
so that X(ip, g,0) = Z(ip,g). □ 

A cost of this representation is that we have to repeat the stability analysis 
for the interpolated interaction. In the linear approximation, this follows 
from an explicit calculation. The following bound is uniform in the interpo- 
lation parameter. 

Lemma 4.6 For all ((f>,g,t) SlxGx [0, 1], we have that 

/d M( i_ t)7 (£)n0 + > |<A 4 -15t; 2 <7. (78) 

Proof 

dM ( i_t) 7 (0 V(<p + tg) = g J dM ( i_t) 7 (C) P^ + i,v) 

g P^v-il-t)^ 



> g iZ--l5[ v -(l-t)<y] } (79) 



and 



1_J_ = v _ 1 < v _(i_ t ) 7 < v = (80) 
show the assertion, since f3 2 = L 2 ~ D < 1. □ 
4.2.1 Large field domination 

We are now in the position to estimate the downstairs factor in ([731), using 
up a fraction, say one half, of the large field behavior of (78). 



Lemma 4.7 Let V be given by (6C). For all ((f), g) £lxG, we have that 



2 

< 



C(g)^g exp(-|</> 4 + 15<?*/ 



(81) 



with 

C(g) < 87 I A 6 + 9v*A 2 g) (82) 
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where 



...4 



A 2n = sup I e""r </> 2n ) . (83) 



Proof The stability bound ([78]), in conjunction with the elementary esti- 
mates 

/ o \ 2 

(— P 4 (0,«)J = 16 (/ -6t>0 4 + 9t; 2 2 ) < 16 (0 6 + 9u 2 2 ) 

(84) 



and 

XP ("4- J " -I 4 

(85) 

implies that 

exp^ dp (1 . ()T (0^ + ^)j 2 (^/ d ^h(£) y ^ + ^)) 



< exp ( -^- + 15vg) 87^ ( D + 9*/ 

< exp ( -i^- + i5 vg \ 8jg? (A 6 + 9gv 2 A 2 ). □ (86) 



Eq. ( pq ) shows that each </> in the downstairs factor kills g^ . Therefore, each 
of the two (^-derivatives in its calculation yields g^. Two (^-derivatives give 
a total factor of 32. 

4.2.2 Fluctuation integral 

Half of the stability estimate has now been used up for the control of the 
downstairs factor. The other half suffices to do the fluctuation integral. 

Lemma 4.8 Let b and c be positive constants. For all ip G R, we have that 
/d.fOe-^^^fl). (87, 

Proof From (^), we deduce that 

^ 4 > b i^- b 4- ^ 

The Gaussian convolution of a Gauss function is again a Gauss function. 
From its explicit form, we find (|87|). □ 
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Lemma 4.9 For all (</>, g) £ R x G, we have that 

d^ 7 (C) e-J"^a-*)-r«) v(^K+€) 2 dAt (a _ t)7 (0 K(0 + C + O 

C( 5 ) 75 exp (- i T --^- 2 + 15 v 2 g + i ^) . (89) 



2 

< 



Proof From (||) and (§3), with b = f and c = 4<" 1 , it follows that 

J ^(Oe-S^-^n^+i)! (J d/i (1 _ t)7 (0 W + C + O 

< y d/i t7 (C)C(g) exp(-|(0 + C) 4 + 15i; 2 <7 



^C^^expl-i^^-^ + lS^g + i^- 1 ). □ (90) 



In eq. (|73|), we also encounter another fluctuation integrals without down- 
stairs factors. It is estimated in the same manner. 

Lemma 4.10 For all {4>,g) G R x G, we have that 

f d^ 7 (C) exp (- J ^ {1 _ t)l {i) V(<j> + ( + 



< 



exp 



( - - ^ + 15v 2 g+- g 2 ?- 1 ) . (91) 

\ 2 1+ gP~i Y y 8 y ' 



Proof From (||) and @, with 6 = f and c = 2g p ~ 1 , we find that 

y d^t 7 (C) expf-yd/i ( i_ t)7 (0 V(0 + C + O 
< y dMt 7 (C) exp(-|</> 4 + 15u 2 5 



* ex K"^rT^^ 2 + 15u25+ ^ 2 " 1 )- D (92) 

4.2.3 Scale transformation 

The remaining task is to combine (|89|) with (|9l|) and to rescale the field and 
the coupling. We insert these estimates into fl73|) to obtain the following 
error bound. 
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Lemma 4.11 Let F\ be given by (68). For all e £ (0, |) and aZZ L £ 

{2,3,4, . . . } there exists a maximal coupling g max > depending on L, e, such 
that for all (ip,g,t) GRxGx [0,1], we have that 



< 



g'-'Ftig) Z QU (<t>,g). 



(93) 



Proof Insert (§g) and @ into @ to conclude that 



^X(^,g,t) 



< a(6g)$ C(5g) exp ( - ^ - (/3^) 2 



a (M p 
2 l + (<5^ 

exp ( 15n r-^fy+ { ^1 + I) (J^) 2 ^ 1 ) . (94) 



We choose <? max such that 



a 5? g e C(5g) exp {15 (a 5 - 1) v 2 g) 



and 



a ^rfCF " bQu{9) - ° (96) 

The condition ( p6[ ) is easy to fulfill because a f3 2 = L 2 is on our side. The 
condition @) is also easy to fulfill, but it requires <7 max to be exponentially 
small as a function on L. 



Lemma 4.12 Let e, L, and g mSbK be as in Lemma Put 

F 2 (g) = g\~* Fxig). (97) 

Then \\T X {Z)\\ F2 < 1. 

Lemma ( [4. 1 1|) is the first instant in this section where we need a small 
coupling argument. (Additionally, this F\ limits the contraction mapping 
to small couplings.) To deal with large couplings case, we have to look for a 
modification of (|32|). Since also (p8[) would require a modification, and since 



the exponent a = % ~ e 1S anyway too small to meet the condition fl47), 



we will not eleborate on this possibility here. Instead, we will replace ( p0| ) 
by an approximant from higher order perturbation theory and modify the 
estimate of this section for this case. 
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5 Perturbation theory in g and g 2 ln(g) 



In this section, we first recall the formal power series solution to the fixed 
point problem T(V) = V, where T is defined by S(Z) = exp(— T(V)) 
with Z = exp(— V). As in |4l|] , we develop V into a double perturbation 
expansion in both g and g 2 ln(g). We then prove a stability bound for the 
perturbative approximants (of odd order in g), extending the analysis in 



a 



5.1 Formal power series representation 

In three dimensions, the renormalized ^-trajectory is not expandable into a 
formal power series in g. However, it does admit a formal power series repre- 
sentation in both g and <? 2 ln(<7). See [^TJ. To simplify the bookkeeping, we 
prefer g and k = ln(<?) (instead of g 2 ln(g)) as formal expansion parameters. 
Let V(4>, g, k) be given by a double formal power series 

oo [2] 

V(M = EE V^M) g r k" (98) 



r=l a=0 



with polynomial coefficients 



N(r,a) 

V^ a \<j>) = E P2n(^v)V^ a \ N(r,a) = r-2o+l. (99) 

n=0 

The maximal number of fields at a given order (r, a) is peculiar to cfi 4 - 
theory. For safety reasons, we define 

v£' a) = n > N(r,a). (100) 

Also, we set the order zero to zero. To first order, the trajectory is defined 
to be a pure normal ordered cj) -vertex, 

v£ a) = Ko s n , 2 . (101) 

The perturbative fixed point turns out to have two free parameters one for 
each resonance. See [fllf . All of these solutions are suitable approximants 
for the contraction mapping. We set both parameters to zero, 

v^ m = y (3 ' 0) = 0. (102) 

The choice ( |102| ) has the advantage is to have a minimal number of vertices. 



The formal power series (98) supplies us with a sequence of polynomial 
approximants 

'"max [ 2 ] 

F (w) (0)5) = £ v {r > a) {4>) g r Hg) a (103) 

r=l a=0 



24 



labeled by the maximal power r max of g. The first of which, r max = 1, is 
the above linear approximant. The default value of rmax will be seven in the 
following. 



5.1.1 Recursion relation 

To be a fixed point of the extended renormalization group transformation 
T (the transformation for the interaction V) in the sense of a double formal 
power series, the coefficients have to satisfy the following recursion relation. 
Let (0\;... ;On) 7< j, denote the cumulants associated with the Gaussian 
moments 

(o 1 ---o n )^ = y"d^(c)0i(* + C)---0n(* + O- (!0 4 ) 



Lemma 5.1 Let V be given by (98). Let KiV)^^ be the coefficients defined 



(-i) 



i+l 



i=2 



EE""EE 5r ^ 

ri=lai=0 Ti=l cii=0 



(2n)lv r > 



dfi v (4>) P 2 n(^,v) (V(^);... ;V ^ a A T . (105) 

\ / <y,pi/j 



Then Z = e v is a fixed point of S in the sense of a formal double power 
series if and only if 



[\-L 



3— n— r\ 



V, 



(r,a) 



2n 



-LT 



E 

b=a+l 



/A HL) b -«v£ b) + K(v)t ] - 



(106) 



To derive this set of equations, one performs a cumulant expansion for the hi- 
erarchical renormalization group, rescales the coupling, and compares equal 
double orders (r, a). 



Lemma 5.2 The system of equations (106) has a unique solution of the 
form (98) with the properties (101) and (10k). 

Proof The set of equations ( |106| ) can be solved recursively. The condition 
( |100| ) iterates through the recursion.]] One proceeds forwards in the order 
r — 1 — > r and, at the order r, backwards in a — * a — 1. Suppose that we 
have computed both 



6 The Gaussian integral projects onto the P2n('i/')-coniponent of the cumulant. 
7 One cannot build connected diagrams with more than 2(i — 2 a + 1) external legs from 
r — 2a vertices g : cj> 4 -v and a vertices g 2 ln(<j) : <fP : v . 
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1. for all (s, b) with 1 < s < r — 1 and < b < [|1 and 

2. for all (s, b) with s = r and a + 1 < b < |. 

Then this data determines the right hand side of 
compute V^n"^ for all n < N(r,a). We find two cases. 

1. Non-resonant case: If 3 — n — r ^ 0, then ( |106| ) determines V^'^ ■ 

2. Resonant case: If (r, n) 6 {(2, 1), (3, 0)}, then the left hand side of 
( |106j ) is zero. In both cases, we find a constraint on the right hand 
side of (|0|). 

The two resonances are resolved by logarithmic corrections. Consider the 
mass resonance (2, 1). Since^] 

K{V)f> 1] = 0, (107) 

the equation labeled by (r,a,n) = (2,1,1) is automatically satisfied. The 
equation with (r, a, n) = (2,0,1) becomes 

= — L~ 2 ln(L) V 2 (2 ' 1} +K{V)f' 0) . (108) 

(2 1) (2 1) 

We use it to determine Vn ■ The other parameter Vq~ ' is unconstrained. 
We put it to zero. The vacuum resonance (3, 0) is analogously resolved. □ 



106). Therefrom, we 



5.2 Stability bound 

To any finite order, perturbation theory furnishes approximate solutions to 
our fixed point problem, which are polynomials in <j) with coefficients that 
are polynomials in both g and g 2 ln(g). We intend to a polynomial of this 
kind as the approximate fixed point in our contraction mapping. To this 
aim, we need to prove two properties, a stability bound and an error bound. 
Both bounds will be proved analogously to [44]. 



5.2.1 Tree approximation 

We first prove stability for the tree approximation, which is defined as the 
polynomial in </>, whose coefficients are simplified to their leading powers in 
g. This bound extends by continuity to a stability bound for the complete 
perturbative approximant in a small coupling region. 

(r) (r 0) 

The set of coefficients V tr J e = V^ir+i) can ^ e com P u ted independently 
of the others. They define a tree approximation 

oo 

v tree {^g) = E ^ 2(r+1) 9 r v};l (109) 

r=l 

8 g 2 ln(g) : <f> 2 : v is not generated in the contraction of two vertices g : cj> 4 : v . 
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to the renormalized </> -trajectory^]. Notice that we have replaced : 4> 2n : v by 
its highest term cp 2n . 

(r) 

Lemma 5.3 The recursion relation for V tr J e decouples. It reads 

1 - L^A V& = 



T 3-r \ ^ A T/ (ri) \A n) 

L 2.^1 1\ / , / , °r,ri+-+n Vfree V tree 

i=2 n=i ri=i 



X (2(r + l) 1 )!^'+ 1 ) f dfhj ^ P %+1)^>^) ^2(r 1+ l);--- ;^2(r i+ l) 



(110) 



All vertices in the tree approximation are irrelevant in the extended power- 
counting. In particular, there are no resonances in the tree approximation. 
The tree coefficients have a simple sign pattern. 



Lemma 5.4 For all r > 1, V^} e = (-l) r+1 



Proof An induction on the order r. □ 



^free 



It follows that all tree approximants with even maximal order r max are un- 
stable. Therefore, we will restrict our attention to the friendly approximants 
with odd maximal order. 

The sign pattern remains valid at sufficiently small couplings. To realize 
this, assemble the perturbative approximant with loop contributions. 

Lemma 5.5 Let V^ree = ^tree = 0. Then we have that 



i max 



+i 



EE^ (r,a W Ms)- = E ^ 2n 9 n - 1 {vlz;^ + hn(g) 

r=l a=0 n=0 ^ 

(111) 

with \2n(g) = 0(g,g 2 ln(g)). 

The tree coefficients are the leading terms of the perturbative vertex func- 
tions at small couplings. The loop corrections are continuous functions of 
g (since they are polynomials in g and g 2 \n(g)). Therefore, properties like 
the sign pattern at g = extend to a finite region g G [0, g max \ of small 
couplings. 

9 The sum of tree graph contributions is in fact convergent. We will not use this fact, 
since we are dealing with finite order approximants, which are polynomials in (j>. 
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5.2.2 Effective 4 -coupling 

All g 2 ln(g)-terms are subleading. For this reason, the tree graph bound of 
[44 1 applies also to the three dimensional model. 

Let ^2n{g) = Vtree ^ + ^2n(g) (a polynomial in g and g 2 \n.(g)) so that 
the perturbative approximant becomes 

' max +1 

V^\<f>,g) = <? n 9 n ~ X ^n{9). (H2) 

n=0 



Lemma 5.6 For all r max > 1, there exists a maximal coupling <? max > 
( depending on , r maXy J such that fov all g £ & and ti £ {2, 3, . . . , r m ax H~ 1}? 

M2nG?) = ("I)" \V2n{9)\- (H3) 

In the following, we will assume that g max is sufficiently small such that 
([ill ) holds. 

The following statements are presumably true for any finite order r max £ 
2N + 1. As a part of their proofs, we will have to compute certain coefficients 
recursively. I have only done this up to the (already ridiculously high) order 

^max — 99. 

Lemma 5.7 Let r max £ {1, 3, 5, . . . , 99}. Then there exists a maximal cou- 
pling <7 max > such that, for all g £ G, (Hi ) is bounded from below by 

i 

V^(<p,g) > Y, ( t >2n 9 n ' 1 ^n(g) + ^g P A(g), (114) 

n=0 

where P4,(g) is the solution to the recursion relation 

P4n(g) = fi4n(g) ~ ~. 7 C (H5) 

4p4 n+i (g) 

with the initial condition 

P2(r max +l)(f) = /*2(tw+i)C»)- ( 116 ) 

Proof The proof is an induction on powers of</> 4 . (Notice that 2(r max + l) £ 
4N.) The induction step follows from 

<t> 4n g 2n ~ l j/Mn(sO + <t> 2 g Pin+2{g) + <t> 4 g 2 pin+i{g) 

±An 2n-l J / \ Pin+2(g) 2 , s ( ,2 . PAn+2 X " 

= 4> g s Pin{g) - -. rx + P-in+dg) [4> g + 



4p 4n+4 (5) \ 2p 4n 



i+4 



(117) 
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since p^n+^g) is positive for small couplings. □ 

A proof of the positivity of the effective 4n -couplings is given below. The 
solution of the recursion relation (11151) is a rational function 



P(g,g 2 Hg)) 
Q{g,g 2 ln(g)) 

where P and Q are polynomials in g and g 2 ln(^). 

Lemma 5.8 Let r max G {1,3,5,... ,99}. There exists a positive number 
c > and a maximal coupling g ma x > such that for all g G G, we have 
that 

Pi (g) > c. (119) 
Proof The value /04(C)) is determined as the solution of the recursion relation 

(2n) 

P4n(0) = VST 1 * ~ , (120) 

4p 4n+4 (0) 

with the initial condition P2(r max +i)(0) = . Therefore, it depends 

only on the tree graph coefficients. An explicit computation shows that 
/94(C)) is a positive number. Since (|118| ) is a continuous function of g, the 
assertion follows. □ 

The effective (/> 4n -couplings at order seven will be listed below. In particular, 
the value of /04(C)) at order seven is 

/ \ 4306 
«(0) = m . (121) 

As a side remark, we mention that the effective t/> 4 -coupling in the tree ap- 
proximation is not a small number at large orders. (It presumably converges 
as the order is taken to infinity.) 

Lemma 5.9 Let r max G {1,3,5,... ,99}. There exist positive numbers 
fi'max; c, and a (all strictly larger than zero) such that for all (4>, g) £lxG, 
we have that 

V^\<p,g) > g (l^-a). (122) 



Proof From ( 114 ) and ( |116| ), it follows that there exist two polynomials A 
and B and a positive number c such that 



V^i^g) > g \A(g,g 2 ln(g)) + B (g , g 2 ln(g)) cj) 2 + c c/) 4 



(123) 
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for all € M and g £ [0 

jflmax]) where G m ax is a certain positive number. 
For this maximal coupling, define ||^4||oo = su P 9 gg \A(g,g 2 ln(^) | and analo- 
gously 1 1 B | \ tx> . Then we have that 



V 



('"max) / 



,g) > 9 

> 9 



\A\\ 



\BWoo6 2 + c6 4 



2c 2 



(124) 



Put a 



IIS II - i 

^4 oo + to obtain the assertion. □ 



2c 



The remaining stability analysis is completely analogous to the linear case. 



Lemma 5.10 Let r max E {1,3,5,... ,99}. Let <? max 6e as in Lemma 5.L 
Let Fx : G — * M + be given by 



Fx(g) = exp ( ag + ^g 



2p-l 



(125) 



Then Z( rmax ) = exp{— y( r max)J j s b oun ded in the norm associated with F\. 
We have that 



I ^(Vmax) I 



< 1. 



(126) 



This shows that the perturbative approximants are indeed in the domain of 
the extended renormalization group. The stability bound is complete aside 
of a proof of the positivity of p^n+i (g) ■ By continuity, it suffices to prove 
the positivity of p4 n+ 4(0), which depends only on the tree coefficients. 



5.2.3 Computation of tree coefficients 

There is another way to compute the tree coefficients than by the recursion 
relation ( |110|) , which uses a Hamilton- Jacobi differential equation. This 
other way is both simpler than to iterate ( |110| ) and it also relies on an 
interpolation formula, which we will need independently in the error bound. 

The perturbative renormalization group is the formal power series solu- 
tion of the non-linear transformation 



T(y)(ij>,g) 



-a In 



d/i 7 (C) expi-VW + CSg)} 



(127) 



It can be computed in two steps. Step one is the fluctuation integral 



In 



d Mf7 (C) exp{-V(<t> + C,6g)} 



(128) 
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evaluated at t = 1. The interpolated quantity satisfies the renormalization 
group differential equation 



o 2 



■W(<f>,g,t) 







W(</>,g,t) 



(129) 



in the sense of a formal powerseries, with the inital condition 

W^,g,0) = V(c/>,g). (130) 

Step two is the scale transformation of the result of step one, 

T(V)^,g) = aW(M,5g,l). (131) 

Consider the tree approximation hereof. The tree approximation affects only 
step one. Eq. ( |12g| ) has to be replaced by the Hamilton-Jacobi equation 



d 

QjW t ree((/>,g,t) 



1 







5 { —Wtree(<P,g,t) 



with the initial condition 

W t ree((p,g,0) = Vtree(<f>, 9)- 

The condition of renormalization invariance becomes 

V t ree(4>,g) = a W tr ee{Plp,5g,i) . 



(132) 

(133) 
(134) 



Lemma 5.11 The Hamilton-Jacobi equation (13i) has a unique formal 
power series solution 



W{<f>,g,t) = J^B 2n (t)<t> 2n g n - 1 



n=2 



with the boundary condition 

W{cp,g,l) = a- l W((3- l 4>,5- l g) 
and -E>4(0) = 1. It reads 



B 2n (t) = b 2n i ~7 



1 



+ t 



n-2 



L 2 - 1 

where the coefficients b 2n are recursively determined by 
(n - 2) b 2n = 2 ^2 ml b 2m hi 

m+l=n+l 



(135) 



(136) 



(137) 



(138) 



with the initial condition bi = 1. 
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We remark that this solution makes sense beyond a formal power series. To 
see this, one writes 

W(<j>,g,t) = g^W^^t) (139) 

and shows inductively a bound on the positive coefficients b 2n - However, 
since we only need the formal power series solution, we leave this issue 
aside. 

Corollary 5.1 The tree graph coefficients are given by 

v£l = b 2(r+1] {j^Jjl) ■ ( 14 °) 

This confirms the sign pattern of the tree coefficients. 

Once the tree coefficients are written on the blackboard, we proceed to 
compute the effective </> 4ra -couplings in the tree approximation. With the 
t-dependence switched on, their recursion relation reads 

P4»(t) = B 4n (t) ~ ( 141 ) 

starting at 

P2(r max +i)(*) = £ 2(rm « +1) (t). (142) 



Lemma 5.12 The effective (j) -couplings are given by 



P4n(t) = r 4n <^ 7 



1 



+ t 



2(n-l) 



with (t -independent) numbers r^ n determined recursively by 

J 4n+2 



r 4n — 



rin+4 



where 



(143) 



(144) 



(145) 



r 2(r max +l) - ft 2(r max +l)- 

Remarkably, the effective (f> -coupling comes out to be independent of the 
interpolation parameter t. To seventh order we find the following numbers. 



Tree coefficients 


n 


(-l) n &2„ 


2 


1 


3 


-8 


4 


96 


5 


-1408 


6 


23296 


7 


-417792 


8 


7938048 



Seventh order 


n 




4 


7938048 


3 


33817/19 


2 


90032/1321 


1 


4306/5627 
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The error bound is complete for the order seven approximant. As a sidedish, 
we find the following useful bound which is uniform in the interpolation 
parameter. 

Lemma 5.13 The tree approximant of order seven 

8 



W^,g,t) = Y. B *n(t) <P 2n 9 n - 1 (146) 

n=2 



satisfies for all ((f), g,t) £ R x R + x [0, 1] the lower bound 



wgL&, g ,t) > (147) 



The recursion relations ( |138| ) and (144) can be solved by computer algebra. 
Their solution proves the positivity of the tree approximation (at least) up 
to the order 99. 



5.3 Error bound 

To prove an error bound for the higher order approximants, we proceed 
analogous to the linear case. The main tool is a generalization of the inter- 
polation formula (76). From perturbation theory, we have a polynomial 



= ^^(^W/ Hg) a 



r=l a=0 

which satisfies the scaling relation 

n 



n=l 

Here 7 ? ^ rmax ^ denotes a projector 
P (rmax) (<f ln( 9 ) a ) = 



T 



c/' r ln(5f) a if r < r max and 
else. 



(148) 



(149) 



(150) 



The truncated cumulant expansion in ( |149| ) contains terms of higher order 
than g rmax . These are projected out by means of , p( rmax ). 

Definition 5.1 Let : R X G X [0, 1] be defined by 



W {r ^(ij,g,t) = - ^ (rmax) ( [^ (nnax) (-^5);] ") 

n=l 



n+1 



n\ T 



(151) 
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This interpolation is identical with the formal power series solution of ( |128| ), 
projected to , p( rmax ) . Its boundary values are 

W( r ™°*\ip,g,Q) = V^ r ^\tp,g) = aW^^^^g,!). (152) 

We use it to define the following generalization of (|70|) (the case r max = 1). 

Definition 5.2 Let X^'^ :lxGx[0,l] be defined by 

X^\ip,g,t) = Ijd^iC) expf-W^ )(J3fj, + C,5g,l-t)\\ 

(153) 

To be well defined, eq. (153) calls for a stability bound for the interpolation 
(|152| ). Postpone this issue for a short while. Eq. ( |153[) yields the following 
representation for the error term. 



Lemma 5.14 Let r max £ {1,3,5,... ,99}. Then ( \15S\ ) is well defined ft 
the perturbative approximant We have that 



or 



^(^^(V'.s) = X (rmax) (^,5,l) -X (rmax) (V>,3,0). (154) 

The usefulness of this representation relies on the following differential for- 
mula. 



Lemma 5.15 Let r max E {1,3,5,... ,99}. Then ( 155 ) is continuously dif- 
ferentiable in t G (0, 1). We have that 

d 



-X^\nl>,g,t) 



a 



dyu t7 (C) expf -W< r ™* )(Jty + C, Sg, 1 - t) 
x f dMfry(C) expf-W^^ + C^l-*) 



a-1 



(155) 



Proof 

di 



d// t7 (C) exp 

d/%(C) 
d/i t7 (C) exp 



7> |ln y d/i (1 _ t)7 (0 exp(-y(^ + C + e,9)) 



5 7 # 



-V 



d 7 3 2 
di + 2dC 



V 



-V 



1 ( 9 
+ - ^zV 



2 \d( 



(156) 
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and 



\dt 2 d( 2 ) \ / 2 dc 2 

= vl (± 

= Hl p { }) 2 - D < 157 > 

The upstairs factor is understood as a synonym for the truncated cumulant 
expansion. It is a polynomial expression. The downstairs factor in ( |156| ) 
is in particular a polynomial expression and (157) is a manipulation of a 
polynomial expression. 

The important feature of ( |155| ) is that in the downstairs factor all orders 
lower than r max are cancelled by the i-dependent upstairs factor. To be 
well defined for all values of the interpolation parameter, ( |153j ) requires an 
additional stability bound. 

Lemma 5.16 Let r max € {1,3,5,... ,99}. There exist positive numbers 
5max, c, and a (all dependent on r max/ ) such that the following stability bound 
holds for all ((f), g, t) G R X G X [0, 1]: 

W^\4>,g,t) > g (l^-a)- (158) 

Proof 

Since W^™^ ((ft, g, 0) = V^ rni ^\(p, g), we know that the bound (|158f ) is 
valid at t = 0. Furthermore, we have shown that the effective (/> 4 -coupling 
is independent of t in the tree approximation. The assertion now follows 
from the uniform continuity of the effective (^-coupling of the complete 
perturbative approximants. 

As the result of a truncated cumulant expansion, we have that 

''"max +i 

\cf>,g,t) = £ ^g n -^ 2n (g,g\t). (159) 

n=0 

For n > 2, each coupling fi2 n {g,g 2 ,t) is the sum of a tree term and loop 
contributions 

V2n{g,g 2 ,t) = B 2n (t) + \ 2n {g,g 2 Hg),t). (160) 



The tree term is given by (137). The loop contributions are higher order 
corrections 

X 2n (g,g 2 H9),t) = 0{g,g 2 \n{g)) (161) 
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as they are polynomials in g and g 2 ln(g) whose coefficients are polynomials 
in t. 

Consider the effective ^-coupling Pi(g,t) defined as above J 10 ! As it is a 
continued fraction of couplings ( |160|) , it is a rational function (of g, g 2 ln(g), 
and t) on some rectangle [0, e^aJ x [0, 1]. Since the tree approximation has 
this particular t-dependence, we have that /54(C), i) = r± for all t G [0, 1] at 
g = 0. Let c = r±/2. By continuity, there exists a positive number <J max 
(with < g ma , x < 5max) sucn that for all (g,t) G [0,g max ] X [0, 1], we have 
that 

P 4g,t)>c. (162) 

Taking care of the constant and quadratic term in (f> analogously to ( |123| ) 
and ( |124| ), the assertion follows. □ 

As in the linear case, the stability bound on the interpolated interaction 
is independent of the interpolation parameter. The remaining analysis is 
completely analogous to the linear case. 

5.3.1 Large field domination 

The harvest of the higher order perturbation theory is a higher power of g in 
the bound after dominating the large fields by part of the stability estimate. 



Lemma 5.17 Let r max G {1,3,5,... ,99}. For all (<j), g, t) G R x G x [0,1], 
we have that 

2 



W {rmax) ^ g,t) \ - 1 - 7>( r — ) ( £7 W^™^ ((/), g, t) 







{-^99 +ag 



exp 

< C(g) g r ^l 2 

for some polynomial C G M. + [g,g 2 ln(g)] (with positive coefficients). 
Proof 

The downstairs factor is a polynomial of the form 



(163) 



W {rmax) (cf),g,t) 



2 L J V<90 

^max~!~l 



n=0 



+ £ 



i2n „n— 1 



g n - V B 2n (g,g z ln( 5 ),t) 



(164) 



10 To be precise, we should consider the collections of all 4> 4n - couplings p4 n (g,t) and 
repeat the following reasoning for all of them. 
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with certain polynomials B 2n . Notice that the projector affected only the 
first sum in (|164j) . Notice also that the highest power of (j> is 2(2(r max +l))— 2, 
where —2 comes from the two (^-derivatives. With the help of the stability 
bound, we find the upper bound 







E 9 rm ^ +1 A 2n \B 2n (g,g 2 \n(g),t)\ 

n=0 

+ Yl 9^" 1 A 2n \B 2n (g,g 2 ln(g),t)\\ (165) 



71 — '/"maxH~2 



where 



A 2n = sup 



ex P I 



kin 



(166) 



Expand the polynomials B and take the supremum of t S [0,1] in each term 
to arrive at the bound ( |163j ). □ 



5.3.2 Fluctuation integral and scale transformation 

The fluctuation integral and the scale transformation are now identical to the 
linear case aside of a minimal cosmetic modification to include the constant c. 
We therefore do not repeat them here and jump to the following conclusion. 



Lemma 5.18 Let r max G {1,3,5,... ,99}. Let F\ be the function j\125\ ) 
from the stability bound. For all e 6 (0, ^rp) and all L £ {2, 3, 4, . . . }, there 
exists a maximal coupling g mSuX such that for (ip,g,t) £ 1R x G x [0, 1], we 
have that 



X^ Tm ^{tp,g,t) 



''max 

< 9 2 



F\{g) Z QU (ip,g). 



(167) 



The error bound is an immediate consequence hereof. 



Corollary 5.2 Let e, L, and g m!lx be as in Lemma 5. It. Put 

F 2 (g) = g T ^- e F 1 {g). (168) 
Then we have that Z^'™^ = exp(— y( r ™)^ satisfies the bound 

|| Tl (Z'W)II < 1. (169) 
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We have computed y( rmax )(g) as a polynomial approximant of a formal 
double expansion in g and g 2 fn(g). Then we have shown that, for sufficiently 
small (but finite) couplings, Z( T ' max ) = exp(— y( rmax )) satisfies both 

1. the stability bound ||Z( rmax ) < 1, where F\ is a function of the form 
(||), and 

2. the error bound ||Ti(Zi)||_f 2 < 1, where is given by (|46|), with 
exponent a = ^f^- — e. 

For r max > 7 and e not too large, all assumptions of the contraction mapping 
are satisfied. The construction is complete. 



6 Conclusions and outlook 

The iteration of the contraction mapping provides a convergent representa- 
tion for the ^-trajectory. It can be used to study the properties of the fixed 
point g). One important problem, which can be shown, but which 

will not be shown here, is that g) is positive. A brief discussion of its 



positivity is contained in [44|. Other questions about Z+(<j),g) could also be 
studied in principal, for instance the summability of perturbation theory, 
and analyticity properties of its Borel transform. 

A very interesting question is the behavior of Z*((j), g) at large couplings. 
Conceivably, Z+(cf), g) connects the trivial fixed point at g = with the non- 
trivial infrared fixed point at g = oo. The contraction mapping is potentially 
capable of a construction, which is uniform in the running coupling. But 
such an enterprise requires a better approximate fixed point Z\{(p,g) than 
the one from perturbation theory. It is conceivable that one could extend 



the approximants from [[31], 32, [33], 34] to achieve this aim. 



The underlying scheme of this paper is to compute a renormalized tra- 
jectory as a renormalization invariant curve in the unstable manifold of a 
renormalization group fixed point. This scheme is certainly translatable to 
virtually every theory treated so far with the renormalization group. In par- 
ticular, all hierarchical models mentioned in the introduction can be handled 
that way. We hope to present an extension of this method to the framework 
of polymer expansions and full models in future work. Another aspect of 
this theory is the question how traditionally computed renormalized actions 
converge to the renormalized trajectory. In other words, what is the domain 
of attraction of this extended fixed point of an extended renormalization 
group. This question is related to the problem of renormalization group im- 
proved actions and also to the question how to truncate a renormalization 
group such as to maintain control of the errors. We hope to make progress 
on these and other questions in this context in future work. 
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